Ab initio derivation of multi-orbital extended Hubbard model for molecular crystals 
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From configuration interaction (CI) ab initio calculations, we derive an effective two-orbital extended Hub- 
bard model based on the gerade (g) and ungerade (u) molecular orbitals (MOs) of the charge-transfer molecular 
conductor (TTM-TTP)l3 and the single-component molecular conductor [Au(tmdt)2]. First, by focusing on the 
isolated molecule, we determine the parameters for the model Hamiltonian so as to reproduce the CI Hamil- 
tonian matrix. Next, we extend the analysis to two neighboring molecule pairs in the crystal and we perform 
similar calculations to evaluate the inter-molecular interactions. From the resulting tight-binding parameters, 
we analyze the band structure to confirm that two bands overlap and mix in together, supporting the multi-band 
feature. Furthermore, using a fragment decomposition, we derive the effective model based on the fragment 
MOs and show that the staking TTM-TTP molecules can be described by the zig-zag two-leg ladder with the 
inter-molecular transfer integral being larger than the intra-fragment transfer integral within the molecule. The 
inter-site interactions between the fragments follow a Coulomb law, supporting the fragment decomposition 
strategy. 



I. INTRODUCTION 

In the long history of research on molecular solids, various 
materials have been synthesized and a rich variety of phenom- 
ena have been discovered, e.g., Peierls insulator, Mott insu- 
lator, charge-ordered state, antiferromagnetic state, spin and 
charge density wave state, and superconducting stated For 
the description of almost all these phases, it has been rec- 
ognized that only one frontier orbital, the highest-occupied- 
molecular orbital (HOMO) or lowest-unoccupied-molecular 
orbital (LUMO) plays a crucial roleP^ The tight-binding ap- 
proach, where a molecule is regarded as a single site and 
the various ways of molecular packing are reflected by the 
anisotropy of inter-site transfer integrals and of inter-site 
Coulomb repulsions, has been successful in the systematic un- 
derstanding of the origin of various phases. 

Recently theoretical efforts have been devoted to ac- 
curate derivation of the single-orbital Hubbard Hamilto- 
nian from density-functional theory (DFT; Refs. ISHTOb and 
wavef unction-based calculations UJJH3J magnitude of the 
bare "on-site" Coulomb interaction within the same molecular 
orbital (MO) of the benchmark TTF molecule was evaluated 
as ~ 5.9 eV, which is reduced to 4.7 eV by taking into account 
the intra-molecular screening effects. 5 The on-site Coulomb 
repulsion is further reduced for larger molecules, as expected, 
e.g., in BEDT-TTF molecule 67 where it is estimated ~ 4.2 
eV. It is generally expected that the intra-molecular screen- 
ing effects are not so pronounced for larger molecules due 
to the nature of delocalized molecular orbitals. 5 13 The inter- 
molecular screen ing e ffects have also been analyzed by sev- 
eral approaches pM±o| 

Multi-molecular-orbital properties have attracted recent at- 
tention in the single-component molecular solids EHUD such as 



[Ni(tmdt)2] and [Au(tmdt)2]. Among these single-component 
molecular s olids , the metallic behavior was first confirmed in 
[Ni(tmdt)2]P^^ In [Au(tmdt)2], antiferromagnetic ordering 
was observed with high transition temperature Taf = HO 
KPSEJ] Interestingly, it has been proposed that in [Au(tmdt)2], 
a new type of antiferromagnetic state, called the intra- 
molecular antiferromagnetic (IAF) state, is realized, in which 
two ligands within each molecule have opposite spins P 2 ^ 
From DFT-based ab initio calculations for these compounds, 
several MOs are shown to contribute to the conduction band 
of the crystal sy stems E33I it has been recognized that the 
frontier orbitals for these compounds can be decomposed 
into several moieties, i.e., fragment MOs, and the electronic 
states have been analyzed using the microscopic tight-binding 
Hamiltonian based on the fragment MOsP3 

The charge-transfer molecular compound (TTM-TTP)l3 
(Refs. l3l and l26jj28l) is now considered as another candidate 
of multi-molecular-orbital system! 29 ! 30 ! In this compound, the 
non-magnetic insulating behavior at low temperature has been 
confirmedP^m A charge ordering reflected by the alterna- 
tion of valence of TTM-TTP molecule along the stacking 
direction has been proposed initiallyP^^ The experimental 
analysis based on Raman- scatteringSSEll and x-ray^ measure- 
ments suggested a new type of charge-ordered state, "intra- 
molecular charge ordering (ICO)", which cannot be described 
by the conventional single-orbital approximation. By per- 
forming wavefunction-based ab initio calculations for the iso- 
lated ionic TTM-TTP molecule, 29 we previously revealed that 
this system has a multiconfigurational character. 

In the present paper, we propose a scheme to build 
up an effective two-orbital extended Hubbard model for 
(TTM-TTP)l3 and [Au(tmdt)2], from ab initio multi-reference 
configuration-interaction (MR-CI) calculations. 39 To the best 



of the authors' knowledge, this is the first analysis of parame- 
ter evaluations for multi-MO-based systems. The band struc- 
ture in the so-called normal phase, i.e., no long-range-ordered 
state, is examined in terms of the evaluated tight-binding pa- 
rameters. By noting that the resulting MOs exhibit bonding 
and anti-bonding character between the left and right moi- 
eties, we transform the two-orbital model into the fragment 
MO picture. 

Some of the results have already been presented in Ref . [30| 
where the band structure for (TTM-TTP)I 3 turned out to be 
consistent with the direct band calculation based on the DFT 
results. Our goal in the present paper is to give an explicit 
derivation to extract model parameters in (TTM-TTP)i3 and 
[Au(tmdt)2] compounds. Indeed, one would like to unravel 
the origins of the ICO state in the former and the IAF state in 
the latter. The present paper is organized as follows. In sec- 
tion II, the characteristic feature of MOs is briefly reviewed. 
In section III, we derive the effective two-orbital Hubbard 
Hamiltonian for the isolated (TTM-TTP)+ and [Au(tmdt) 2 ] 
molecules, and evaluate the parameters by using ab initio cal- 
culations. In section IV, the inter-molecular interactions are 
evaluated by focusing on two neighboring molecules in the 
crystal. In section V, we examine the band structure of the 
crystals. Section VI is devoted to the application to the frag- 
ment decomposition of MOs and examinations of its validity. 



II. MOLECULAR ORBITALS 

TTM-TTP stands for 2,5-bis[4,5-bis(methylthio)-l,3- 
dithiol-2-ylidene]-l,3,4,6-tetrathiapentalene and tmdt stands 
for trimethylenetetrathiafulvalenedithiolate. The molecular 
structures are shown in Fig. [T] The crystal structure of 
(TTM-TTP)I 3 is triclinic and the space group is PI: the 
inversion center is located on the mid point of the TTM-TTP 
molecule. The formal charge of the TTM-TTP molecule 
is +1 due to the presence of mono-counterion I3 , and the 
HOMO is singly occupied, that is SOMO (singly-occupied- 
MO) in nature. For [Au(tmdt)2], the crystal structure is 
also triclinic and the space group is PI, where the inversion 
center is located on the metal center. Note that the isolated 
neutral molecule [Au(tmdt)2] is a radical species. All our 
ab initio calculations were performed using the MOLCAS 
7 package. 40 The atomic coor dinate s are read from the 
data of x-ray structure analysis EED Due to the presence 
of inversion center at the mid point of the molecules, the 
resulting MOs can be classified into gerade (g) or ungerade 
(u) with respect to the inversion center. The energy spectrum 
was calculated using a restricted open- shell Hartree-Fock 
(ROHF) procedure. Along these ROHF calculations, three 
electrons are likely to occupy the frontier g and u MOs. 
These MOs are likely to be singly occupied or doubly 
occupied in the ROHF calculations. In order to estimate 
the energies of the expected levels, we set the occupation 
numbers for both MOs as 1.5. 29 Under these conditions, for 
(TTM-TTP)+, we obtain that the SOMO is u and HOMO-1 
is g, in agreement with extended Hiickel calculations,^ and 
the resulting energy difference obtained from the ROHF 
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FIG. 1: The molecular structures and quasi-degenerate molecular 
orbitals for (TTM-TTP) + (left) and [Au(tmdt) 2 ] (right). For (TTM- 
TTP) + , the singly-occupied-molecular-orbital (SOMO) has unger- 
ade (u) symmetry with respect to the inversion center, while the 
second-highest-occupied-molecular-orbital (HOMO— 1), has gerade 
(g) symmetry. For [Au(tmdt)2], the inversion center is located on the 
metal (Au) center, and the SOMOs (HOMO— 1) are of g (u) symme- 
tries. 

calculation is ~ 0.4 eVp^For [Au(tmdt)2], we obtain that the 
SOMO is g and HOMO— 1 is u, in agreement with DFT-based 
ab initio calculations.^ The notations for label of MOs are 
adopted to those in Ref. 22 The energy difference obtained 
from the ROHF calculation for isolated molecule is given by 
(ssomo — £homo-i) — 0.26 eV. Incidentally, the LUMO and 
the third-highest-occupied-molecular-orbital (HOMO— 2) are 
sufficiently separated in energy, (slumo — ^somo) — 6.07 
eV and (e S OMO - ^homo-2) ^ 2.03 eV. 

III. INTRA MOLECULAR INTERACTIONS 

We shall now focus on g and u MOs, shown in Fig. [T] 
and derive the effective two-orbital Hubbard Hamiltonian for 
the isolated TTM-TTP+ and [Au(tmdt) 2 ] molecules. All the 
model parameters are determined so as to reproduce the en- 
ergies of different electronic configurations obtained from the 
ab initio calculations. 



A. Construction of the two-orbital Hubbard Hamiltonian 

In this subsection, we derive the expression of the two- 
orbital tight-binding Hamiltonian for the isolated molecule. 
The relevant orbitals are written as ip g (r) and ip n (r), respec- 
tively. Due to the inversion symmetry around the midpoint 
of the molecule, the wavefunctions obey ip g (r) = ip g (—r) 
and (p u (r) = —p u (—r), where the origin of the coordinate is 
taken as the inversion center (Fig.[T]). 

By focusing on the frontier orbitals, SOMO and HOMO— 1, 
the full Hamiltonian for the isolated molecules is expressed as 

^l-mol = ^1 ^1 £ a C a,(j C a,(j 
a a 

a>,(3,ct f ,(3' a, a' 
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const., 



(3.1) 



where a stands for the creation operator of electron in MO 
a(= g, u) with spin <j(=t,i). The orbital energies for the 
isolated molecules are represented by e^. For the electron- 
electron interactions, we adopt the notation for the integrals 
as 



(aa', /?/?') = J J dndr 2 

x Pa(nW(n)- 



ri - r 2 



■<p*p(r 2 )<Pi3>(r&,2) 



where all the wavef unctions are real. Under the permutation 
of ri and r 2 , we have the relations: 

(aa',/3/3') = (/3/3 1 \aa') = (aa' f3' ' p) 

= (a f a : /3/3 f ) = (a f a,/3 f /3). (3.3) 



Due to the g/u symmetry of the MOs and relations ( |3.3| ), the 
number of independent interactions can be reduced. Depend- 
ing on the choice of the indices a, /?, a' ', /3' in Eq. (3.2), we 
can classify those integrals into intra- and inter-orbital inter- 
actions. The intra-orbital interactions (aa,aa) represents the 
magnitude of interaction for electrons within the same MOs, 
in other words, the "on-site" Coulomb repulsions. Here we 
define the parameters: 



U g = (gg, gg), U u = (uu, uu). 



(3.4) 



For the inter-orbital case, the interactions involving odd num- 
ber of ungerade MO, e.g., (gg, gu) vanishes due to the sym- 
metry constraint. The possible interactions are (aa, /?/?), 
(a/3, /3a), and (a/3, a/3) with a ^ f3, which correspond to 
the Coulomb integral, exchange integral, and pair-hopping in- 
teraction, respectively. We note that from Eq. ( |3.3| ), the am- 
plitudes of the exchange integral and the pair-hopping inter- 
action become identical (a/3, /3a) = (a/3, a/3). We define the 
two-independent coupling constants given by 



J H = 2(gu,ug) = 2(gu,gu), 



U' = (gg,uu) - -(gu,ug), 



(3.5a) 



(3.5b) 



where Jh represents the Hund exchange coupling including 
the pair-hopping term, and U' is the inter-orbital Coulomb re- 
pulsion. 

The level energies are determined by the eigenv alue s of the 
ROHF. We note that the level energies in Eq. (3.1 ) do not 
correspond to the eigenvalues of ROHF, since the Hartree con- 
tributions arising from electron-electron interactions are not 
taken into account. These effects can be incorporated in terms 
of the model parameters U g ,U u , and U' , and we can re-define 
the MO level energies by 

s° g = 4 + £/ g <4 )t c g)t > + V 5>t iCT c UiCr >, (3.6a) 

G 

el = e° u + C/ u (< t c Uit ) + U' $>t, ff c g)(T ), (3.6b) 



where the expectation values are determined by the Hartree- 
Fock calculation. Since the average charges on the each MO 



§/7'and e° u 



were set to 3/2 we obtain s® 

In order to express the Hamiltonian in a compact form, we 
introduce the density operators in the normal-ordered form: 



C g,cr C g,cr 
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4' 



°U,cr°U,cr 



4' 



(3.7) 



andn g = (n g ^ + n g ^),n u = (n u ^ + n u ^). The Coulomb in- 
teractions having amplitudes U g ,U u , and U' can be expressed 
in bilinear forms of these density operators. In terms of model 



parameters defined in Eqs. (3.4), (3.5), and (3.6), and also in 



terms of density operators given in Eq. ( 3.7 ), we obtain the ef- 
fective two-orbital tight-binding Hamiltonian for the isolated 
molecules: 



-Hl-mol 



+ U g n gt n gi + U u n ut n ui + U'n g n n 



(3.8) 



where S a is the spin operator given by S a = 
\ J2 a a i c \x a cr (T,cr / Ca,(T / with ct being the Pauli matrix. 
The energy constant Eq and the level energy e g (e®) corre- 
spond to the total energy and the MO level energy for cp g 
(ifu), obtained by the ROHF calculations, respectively. 



B. Possible configurations 

In this subsection, we introduce the possible configurations 
for the two-orbital system including 0, 1,2, 3, or 4 electrons. 
Let us stress that the MOs used to describe the different states 
remain unchanged and correspond to the three-electron case. 
For this purpose, we introduce the following notation: 



The zero-electron state is given by 



(o)\ _ 



= |o>, 



(3.9) 



(3.10a) 



where |0) represents the Slater determinant in which the MOs 
are occupied up to the third-highest-occupied-MO. 

For the one-electron states, we have two configurations 



t /-4,ti°>< 



= 4 it |o>. 



(3.10b) 



(3.10c) 



Since the same relations hold for the spin-down states, we fo- 
cus only on the spin-up states hereafter. 
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For the two-electron states, there are four possible config- 
urations, which can be classified depending on the symmetry 
and spin states. For the two-electron states with gerade sym- 
metry, we have two configurations: 



D. Parameter mapping 



Id 



g,s,2/ 



n 
n 



(3.10d) 
(3.10e) 



where both configurations are spin singlets. For the two- 
electron states with ungerade symmetry, there are also two 
configurations: 



i 



t C u,f|0)> 



(3.10g) 



where the suffixes s and t represent the spin- singlet and spin- 
triplet states, respectively. The S z = 1 component is only 
considered for the spin- triplet state. 
The three-electron states are given by 



i^ 3) ) 
i^ 3) ) 



u 
t 

t 

u 



4,t c u,t c uJ°>' O.lOh) 

c Lt c u,t c Li°>- 



• f J J ln ' (3.10i) 



Finally the four-electron state has a unique configuration: 



l^ (4) > 



n 
n 



1" 



c g,t c u,T c g,; 



< 4 |0>. (3-lQj) 



C. Ab initio calculations 

The energies of the different configurations were eval- 
uated by performing MR-CI calculations.^ Several ba- 
sis sets were used for convergence control, however, a 
weak dependence on the choice of basis set is observed. 
Throughout this paper, we adopt the basis sets contractions 
for the elements S(7s6pld)/[4s3pld], C(5s5pld)/[3s2pld], 
Au(13sl0p9d6f)/[5s4p4d2f], and H(3s)/[ls]. The SOMO and 
HOMO— 1 levels are well separated from the other MOs. 
Therefore, these two MOs are used to generate a so-called 
model space containing three electrons in two MOs. In order 
to evaluate full parameters, we also consider virtual states by 
removing/adding electrons. 

We note that the two-electron gerade wavefunctions are 
given by the superpositions of the two configurations [Eqs. 
( |3.10d| ) and ( |3.10e| )]. By taking advantage of the information 
conveyed by the wavefunctions, we can access the Hamilto- 
nian matrix expressed in the basis of Eqs. ( |3.10d|) and ( 3. lOe ). 
Detailed formulation is given later (see Sec. IV Q. 



The model parameters in Eq. (3.8) can now be evaluated 
by relating the ab initio calculation results with the respective 
energies expressed in terms of the microscopic model param- 



eters. From the Hamiltonian Eq. (3.8) and the configurations 



[Eq. (3.10)], the energy expectation for the zero-electron state 

E (o) = (^(0)1^.^1^(0)) can be expressed as 

In the similar way, the respective energies 22 W = 
(^(*) \ Hi- mo \\ip^} can be expressed as 



4 d 



E. 



(2) 
g,s,2 



1 1 c 



(2) 
u,t 



E (o 

E (0 
E (0 
E (0 

E (0 

E (0 
E (0 
E (0 



el-\u & -\u>, 
el-\u.-\u>, 



2s° g --U g -3U>, 



2£° - \u u - 3U', 



el+el-\u z -\u^-2U' 



e° + e° 



2s» 



(3.12a) 
(3.12b) 
(3.12c) 
(3.12d) 

(3.12e) 
'-U u - 2U' - \ J H , 



(3.12f) 

\u g - l -U n - 5 -U', (3.12g) 



K + ^-\ U z-\ U "-\ U 'i ( 3 - 12h ) 



2e° g + 2e° u - |t/ g - hj* - 2U'. (3.12i) 



Seven model parameters are to be determined, while ten en- 
ergies are calculated for the respective configurations. Thus, 
the expressions for the model parameters in terms of ab ini- 
tio energies are not expressed uniquely, but we can obtain the 
same numerical values irrespective of the expressions. One 
possible way to relate the configuration energies to the model 
parameters is given by 



En 



■ l jE ;(o) + 3 ^(i) + 3 ^(i) 
16 8 g 8 u 



16 



1 



-E ( 



(i) 



3 (2) _3 ( 2 ) 
3^,8,1 g^g,s,2 



s" = -E^ - IeP 



-25< 4 > 



£ (4) , (3.13a) 



(3.13b) 



2 u 
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-|42i + l4S a + |^ (4) . ( 3 - 13c > 

C/ u = £7<°> - 2£« + 45,2, (3.13e) 

^ = ^ (0) -i45i"^- 2 + ^ (4) ' (3-13f) 
J H = 45-45- (3-13g) 

We note that Jh can also be derived from the off-diagonal 
pair-hopping term 2(^ il \H\^l 2 ). 

E. Results 

The model parameters obtained in this way are summarized 
in Table [I] The on-site Coulomb repulsions U g and U u are 
comparable to the values for the BEDT-TTF molecule ~ 4.2 
eVPZl We have also applied the present scheme to the TTF 
molecule and obtained large value of the "on-site" Coulomb 
interaction ~ 6.2 eV. By comparing the values of U g and 
U u for TTM-TTP+, we find a relatively smaller value for 
U g , which can be explained by noting the left-right bond- 
ing character of the u MO while the g MO is anti-bonding, 
as seen in Fig. [T] On the other hand, the values of U g and 
U u for [Au(tmdt)2] are almost degenerate as a consequence 
of the enhanced derealization onto a larger system. The 
magnitudes of inter-orbital interactions, U' and Jh, are large 
compared with transition metal atom situations, satisfying 
£4tom = ^atom + ^Hatom- In contrast with the transition 
metal atoms displaying a centrosymmetric potential, there is 
no constraint relation among the couplings, U g , U u , U f , and 
J H . Our results show that (U g + U u )/2 « 3 ~ 4 eV and 
(U' + Jh) ~ 6 eV. As a simplest and most intuitive exam- 
ple for the parameter evaluation of two-orbital system, we can 
consider the hydrogen molecule H2 where the full CI calcu- 
lations can be performed and we can get an insight into the 
large Hund coupling in molecular systems.^ 



TABLE I: Estimated parameters for the (TTM-TTP) + and 
[Au(tmdt)2] molecules. All energies are in eV. 





TTM-TTP+ 


[Au(tmdt) 2 ] 


Level ener^ 




-8.63 


-5.40 


Level ener^ 


iys° u 


-8.21 


-5.66 


u g 




3.70 


3.49 


u u 




3.90 


3.45 


u' 




2.82 


2.48 


Jh 




3.19 


3.65 



IV. INTER-MOLECULAR INTERACTIONS 

By extending the analysis of Sec. Ill, we evaluate the 
inter-molecular interactions by focusing on two neighboring 
molecules, i.e., dimer, in the crystal. 

A. Construction of the two-orbital extended Hubbard 
Hamiltonian 

A similar strategy can be applied to obtain the expression of 
the two-orbital Hamiltonian for the two-molecule dimer sys- 
tem. The molecule index is given by I and II. Each molecule 
has gerade and ungerade MOs reflecting the symmetry with 
respect to the each inversion center. For the gerade and unger- 
ade MOs on molecule I, we adopt the notations (p gy i and <^ Uj i, 
respectively (^ g ,n and <^ Uj n for molecule E). The MOs cor- 
respond to the strictly non-interacting situation. Practically, 
they were determined by pulling apart the two molecules. The 
respective MOs are shown in Fig. [2] The inversion transfor- 
mations with respect to the inversion centers located on the 
molecule (Ri and Rr) are represented as, e.g., y? g) i(Hi+r) = 

<£g,i(-Ri — r ) an d <Ai,i(-Ri + r ) = — <Ai,i(-Ri — r )- From the 
x-ray structure analysis, it has been shown that the system ex- 
hibits another inversion center located at midpoint between 
two molecules ! 19 l 26 l For the inversion transformation around 
the midpoint, the wavefunctions of the two molecules follow 
the relations 

^ g ,iO) = <p g ,n(-r)> WiO) = -<Ai,n(-r). (4.1) 

The <£ gj i, y? Uj i, <£ gj n> and <^ Uj n MOs are used to define the 
model space. 

For the inter-molecular interactions, the most important pa- 
rameters are transfer integrals between the molecules. We de- 




FIG. 2: Side view of the "effective molecular orbitals" for the two- 
molecule dimer system. The MOs of each TTM-TTP molecule are 
the same as those calculated for the isolated molecule. The wave- 
functions Lp g i and <£ui (^ g 2 and (^2) are given by the linear combi- 
nations of two g (u) MOs obtained from the isolated-molecule calcu- 
lations. 
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uMO 



gMO 



where a,/3 = g or u, and the molecule index is given by 
hj>i\j'(= Iorll). We assume that all the wavefunctions 



tan 



molecule I 



molecule II 



FIG. 3: Tight-binding parameters for the two-molecule dimer sys- 
tem. We note t gu = — t ug due to the presence of inversion center 
between two molecules. 

fine such transfer integrals (a, /? = g or u) as: 

tap = -(ai\H\/3u) 

= - J dr<p* atI (r)H<p PtlL (r). (4.2) 

In the two-orbital system, four kinds of transfer integral are 
possible (Fig. [3]), which are explicitly given by 



-<gi|#|gn>, 



-(nil^lun), 



-(gil^lun), 



-<ui|fr|gn). 



(4.3a) 



(4.3b) 



(4.3c) 



(4.3d) 



The Hermitian condition indicates (ai\H\/3ji) = (f3ji\H\ai). 
By applying the inversion transformation given by Eq. ( |4.1| ), 
we find t gu = —t ug . 

In addition to the transfer integrals between the molecules, 
we have an extra one-electron interaction term in the two- 
molecule system. If one picks up two molecules from the 
crystal, the inversion center within each molecule is lost and 
the mixtures of MOs, i.e., (gi|if|ui) and (gn|iJ|un) become 
non-zero. Such a symmetry -breaking interaction is given by 

H' A = Ae gu ^( c 1i^ c u,i,<t - 4,n,a c u,n,a + h.c). (4.4) 



In terms of the ab initio Hamiltonian H , this coupling constant 
accounting for this symmetry-breaking interaction would be 
represented as Ae gu = ^(gi\H\\ii) — ^(gn|#|un). In the 
crystal, such interaction disappears due to the periodicity of 
the interactions. 

In analogy with Eq. ( |3.1| ), the electron-electron interaction 
terms are expressed as 



(oiiQii^Pj/Sj,) =11 dridr 2 

(4.5) 



are real and the relation of Eq. (3.3 ) still holds in the present 
case. Under the inversion transformation (r — >> — r) given 
by Eq. ( |4.1[ ), we find several relations, e.g., (giui,uiui) = 
— (gnun, unun). Since the MOs are localized on each 
molecule, we find that the "two-center" integrals given by, 
e.g., (aia{[, /3n/?i) and (ajaj, /3i/3 n ) are negligibly small due 
to small overlaps and cannot be determined within the nu- 
merical error. The interactions expressed as (aia[, 
(a n a4,/3 n /3n), Ai^i), and (^n^n 5 PiPi) are the only 

relevant interactions. We also note here that the hopping pa- 
rameters, Eq. ( |4.2| ), have been defined in the one-electron pic- 
ture and some of two-electron integrals, e.g., (gigi, uiun) and 
(giui, uiun), can contribute to them in the present two-orbital 
systems. However, these contributions are small two-center 
integrals and thus the simple definition of the transfer inte- 
grals can be justified. Depending on the choi ce of indices 
a, a', /?,/?' and z,i', j, j', the interactions (4^5) can be clas- 
sified into five categories: (i) all the indices are identical, (ii) 
three-indices are identical, (iii) two pairs of identical indices, 
(iv) all indices are different, and (v) one pair of identical in- 
dices. We here examine each case separately. 

(i) The interactions are nothing but the local on-site 
Coulomb interactions, given by U g = (gigi,gigi) = 
(gngn,gngn) and U u = (mui,uiui) = (u n u n , u n u n ), which 
are the same relations as Eq. ( |3.4| ). 

(ii) The possible interactions are (gigi, giui), (uiui, uigi), 
and the interactions with the molecular index I and II inter- 
changed. These couplings vanish due to the inversion trans- 
formation within each molecule. 

(iii) In addition to the intra-molecular interactions Jh and 
U', we introduce the couplings: V gg = (gigi, gngn), V uu = 
(uiui,u n u n ), "^gu = (gigi,u n u n ), and Kg = (uiui, gngn), 
where V gu — V ug . These represent the "inter-site" Coulomb 
repulsions. 

(iv) In this category, the only relevant parameter is given by 

I = (giui,gnun) = (giui,u n gn) 

= (uigi,ungn) = (uigi,gnu n ). (4.6) 

This term can be regarded as an "orbital exchange" between 
neighboring molecules. 

(v) In this category, we have only two independent param- 
eters, given by 

X g = (gigi,gnun) = (gigi,ungn) 
= -(gngn, giui) = -(gngn, uigi), 

(4.7a) 

X u = (uiui,g n u n ) = (uiui,u n gn) 

= -(unun,giui) = -(u n u n ,uigi). (4.7b) 

As in the case of the isolated molecule, we introduce 
the density operators in the normal-ordered form: n gi i ja = 

( c l,i,a c g,i,a ~ I) and = (4,i,a c u,i,a - f) where 

i = I, ft. Thus, in terms of model parameters, the full ex- 
tended Hubbard Hamiltonian for the two-molecule system is 
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expressed as 

^2-mol = ^ Hi- mo \ + ^ (A£ g n g) j 

j=i,n j=i,n 



A£ u ra UJ -) 



cr 

- ^gll 5Z( C g,I,cr C U,n,cr ~~ C u,I,cr C g,n,cr + ^-C.) 
cr 

u,i^u,n 

+ V gu (n g ,in u ,n + nu,i%,n) 

( c li^ Cu ' I ^ c g,n,<7' c u,n,a' 
+ Cg,i,<T C u,i^4,n, CT ' c g,n, CT ' + h.c.J 
+ X s J2 [ n g.i( c g,n,<T C ".n >( r + h.c.) - (I H) 

(7 

+ Xu [ n u,i( c Ln,a c u,n,a + h.c.) - (I O E) 

cr 

+ + A£ , (4.8) 
where i^i- mo i is the Hamiltonian for the isolated molecule 



given in Eq. (3.8). The Ae g and As u terms represent the 
energy-level shift due to the neighboring molecules, i.e., 
so-called crystal-field effect. These terms also include the 
Hartree contributions arising from the intermolecular density- 
density interactions, V gg , V uu , and V gu , as in the isolated 
molecule [Eq. ( |3.6| )] . The Hartree contributions from the X g 
and X u interactions can be included into the H' A term [Eq. 
(4.4[)]. The term AEq represents the constant energy shift. 



B. Possible configurations 

For the two-molecule systems, four MOs must be consid- 
ered. Due to the presence of the transfer integrals between 
the molecules, the wavef unctions for the two-molecule sys- 
tem are given by the linear combinations of MOs for the 
isolated molecules. The fragment MOs which exhibit g and 
u characters within each molecule, namely (ty? g ,i, </? u ,i) and 
(^g,n, <Ai,n)» were used to build the symmetry-adapted MOs, 
given by 



(f g 2 



Pu2 

and drawn in Fig. [2] 



V2 
1 

71 
i 

71 
i 

71 



(^g,I + <£g,Il), 

(-<Pg,i + <Pg,n), 
+ ^u,n), 
(<Ai,i + <Ai,n), 



(4.9a) 
(4.9b) 
(4.9c) 
(4.9d) 



As in the case of the isolated molecule system, different 
configurations were considered. In order to specify each con- 
figuration on the new MO basis [Eq. ( |4.9| )], we introduce the 
following: 



g2 U 2 
gl Ui 



(4.10) 



In the two-molecule systems [(TTM-TTP)~]2 and 
[Au(tmdt)2]2, six electrons are likely to occupy these 
four MOs. In order to estimate the magnitudes of the 
Coulomb interactions, we have to add/remove electrons from 
the six-electron ground state. In the present analysis, we 
consider all the configurations with n = 0, 1, and 2 electrons. 
The extensions to the states with n > 3 are straightforward. 
However, we will see that the information is sufficient to 
evaluate the model parameters. 
The zero-electron state is given by 



(o)\ _ 



while the one-electron states are 



= |0>, 



(4.11) 



C T iit |0), (4.12a) 



c L, t l°)' 



(4.12b) 



(4.12c) 



(4.1 2d) 



and c U2 ^ represent the annihilation 



where c gua c Uu(T , c 

operators corresponding to the wavefunctions, cp g i, cp u i, 
and (p U 2, respectively, given in Fig.[2]and Eq. (4.9). 

We also consider all the two-electron configurations. All 
the resulting states are listed in Table [TTJ depending on the g/u 
and singlet/triplet natures. 



C. Ab initio calculations for the two-molecule dimer systems 

The molecule pairs which we consider are extracted from 
the x-ray data, and are shown in Fig. [4] In the case of [001] 
two-TTM-TTP-molecule pair, the "molecular orbitals" ob- 
tained from the ROHF calculation are given in Fig. [2] The 
assumption that we made here is that the MOs obtained in the 
isolated-molecule calculation are similar to the MOs in crys- 
tal. 

MR-CI calculations were performed to extract full parame- 
ters by taking advantage of the wavefunction information. In 
this scheme, one can access the off-diagonal elements of the 
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TABLE II: Possible configurations for the two-electron states in the two-molecule dimer systems, the corresponding expressions in the 
molecular index I and n, and the energy expectation values represented in terms of the parameters in the extended Hubbard Hamiltonian 

El 2) = (^ 2) |#2-moi|^ 2) ), where £ (2) = £ (0) + e g + e u - S(U g + C7 u )/4 - 5U' /2 - 3(V gg + 7 uu )/2 - 5V gu /2. The suffix s and t 
represent the spin-singlet and spin-triplet states, respectively. 



Configurations 



Representation on the molecular-based fragment MOs Energy expectation values 



(2) 



Gerade, singlet '^ g ' s ' 2 



Ida) 



1^,4) 



Ids) 



t4 

4 

t 

n 



n 
t 

44 



\ (I gigi) + I gngn) + I gign) + | gngi}) 

-573d giui) + luigi) - |gnun} - |u n gn)) 

+ 5t7f (Igiun) + |u n gi) - |gnui) - |uig n » 

I (| U1U1) + I untin} - | uiun) - | unui)) 

\ (I gigi} + I gngn) - I gign) - | gngi)) 

~275(l gia ') + I Ulgl ) ~ I g n0lt ) ~ I u ngn» 

-^(|giun> + |u n gi) - | gnui) - | uign)) 

\ (| uiui) + I u n un) + | uiu n ) + | u n ui}) 



£ (2) +£ g -£u-|(7 g + |(7u- 

— V gg + 2^uu — 2^g u _ ^§ 
E (2) + |J H - U+ (tun -tgg) 



£ (2) - £ g + e u + \U S - \U U - \U' 
+f V gg - Ku - l^gu + 2t 

s (2) +£ g -£u-|r/ g + fr/ u - 

-Vgg+lVuu - hV su + 2t 



uu 
17T' 



+ I Jh — \l — (t U u — tgg) 



£ (2) - e g + £u + f f/g - \U» - \U' 



9 ^gu 2t u 



Gerade, triplet 



|(|giui> - |gnu n ) - |giu n ) + | gnui}) 



\{\ gmi) - I gnu n ) + I giu n ) - | gnui}) 



£ (2) " |J H + \I + (tuu - tgg) 



S (2) -|j H + |/-(tuu-tgg) 



Ungerade, singlet 



Ida) 



l^ 2 L) 



t 4 

- 4 
t - 
t - 

- 4 
t 4 



75 (I gigi) - I gngn)) 



£ g -£ u + i[7 g + |J7 u - \U' 



V2 



' 2 V SS 



9 ^gU 



575G giui) + I uigi) + I gntin) + | u n gn)) 

+ 575(1 giun) + |u n gi) + I gnui) + |uig n » 

575 (I giui) + I uigi) 4- I gn un) + | u n gn)) 

-575(1 giun) + |u n gi) + |gnui) + | uign)) 

-75 (|uiui) - I untin)) 



£ (2) + |Jh + |/ -{tuu+t g 



£ (2) + §JH + |/+(iuu+tgg 

E (2) - £ g + e u + §f/ g + iU u 



o V\l\ 



2 v s u 



Ungerade, triplet 



l^L) 



1^,2 



IV$,3> 



l<L) 



t t 

- t 
t - 
t - 

- t 
t t 



£< 2 > + i 



2 K SS ^ 2 



2 v e u 



|(| giui) + I gnun) + I giun) - | uign)) 



|(-| giui) - I gnun) + I giun) - | uign)) 



I uiun) 



-\Jn-\l- (t uu + t g 



£ (2) -§J H -§/ + (tuu + t g 



E {2) -£ g + £u + |t/g 



4 



2 17 



+lv gg ■ 



9 V gu 



Hamiltonian from the knowledge of the eigenvalues space EDGE! 
and wavefunctions. Such approach was originally developed The MR-CI wavefunctions can be written as 
for the situation with larger CI expansions in which effective 

parameters can be extracted by projecting those onto a model = dij\ipi), (4.13) 



9 



(TTM-TTP)I 3 



(a) 



j& jp jr £ 




[Au(tmdt) 2 ] 




FIG. 4: Crystal structure of (TTM-TTP)I 3 [(a) and (b)], and crys- 
tal structure of [Au(tmdt)2] [(c) and (d)]. The target pairs of two 
neighboring molecules are shown by the arrows, e.g., [001] and [002] 
in (TTM-TTP)I 3 represent the nearest-neighboring and next-nearest 
neighboring molecule pairs along the stacking c direction. 



where the |^)s are the configurations defined in Sec. IV B. 
From the MR-CI calculations, the energy eigenvalues {Ej} 
and eigenf unctions {dij} were obtained not only for the 
ground state but all the excited states contained in the model 
space. The Schrodinger equation can be formally expressed 
as 

/ d\j \ ( dij \ 



H 



d 2j 



d 



2j 



(4.14) 



\ d Nj J V d Nj ) 

where j = 1, . . . , N, with N being the number of possible 
states. From we define the unitary (orthogonal) matrix: 

/ d n d 12 ... d 1N \ 

d 2 \ d 22 d 2 N 



U 



(4.15) 



\ djyi dwN J 

Together with the eigenvalues {Ej}, the Hamiltonian can be 



written as 



(e 1 



H = U 



\ 



E 2 



(4.16) 



V E * J 

Thus we can get not only the energy of each configuration, but 
also the off-diagonal matrix elements of the Hamiltonian. 
D. Parameters mapping 

Evidently, the full Hamiltonian is expressed on the basis 
(gi, ui, g2, U2) whereas the model Hamiltonian is written on 



the molecular-based fragment MOs, (gi, ui, gn, un). As can 
be seen from Fig.[2]and Eq. (4^9 ), the correspondence between 
the basis results in 



°gi,cr 



1 

71 

1 
71 

1 
71 

1 
71 



(c g ,i,a + c g>n ,a), (4.17a) 



(-c g ,i,a +c g , n ,a), (4.17b) 
(-c u ,i,a + c u , n ,a), (4.17c) 
(c u ,i,a + c u ,n,a). (4.17d) 



By using these relations, we rewrite all the configurations into 
the molecular-based fragment MO basis, and derive the ex- 
pression of Eq. ( 4.16|) i n terms of the parameters of the model 
Hamiltonian [Eq. ( |4.8| )]. 

The zero-electron state is l^ )) = | 0), and the correspond- 
ing energy can be expressed in terms of model Hamiltonian as 
£7(0) = (^°)|if 2 . mol |^ )) = (2£7 + AE ) - 3(e g + e u ) + 
9(U e + U U + 4U')/8 + 9(V gg + Ku + 2^ gu )/4. 

By using Eq. ( |4.17| i, the one-electron states given in Eq. 
< 4. 1 2|> are re-expressed as 



1 

71 
1 

71 
1 

71 

^(|ui) + |un)), 



(|gi) + |gn», 
(-|gi> + |gn», 
(-|ui) + |un)), 



(4.18a) 
(4.18b) 
(4.18c) 
(4.18d) 



where we have used the notation |gi) = ^|0) and so on. 

By focusing on these four one-electron states, the extended 
Hubbard model Hamiltonian ( |4.8| ) can be expressed by the 
4x4 matrix. The molecular-based fragment MO picture gives 
the following CI matrix expressed in the model space: 



2-mol 



A£ gu 



-As 



As 



gu 



tgU 

M 1} 



gu 



tgu 

-As gu 

^uu 

E (l) 



\ 



Igl) 

|gn) 

|un) 
(4.19) 

3C/ g /4 - 3U'/2 - 3V gg /2 - 

- W u /4 - W/2 - 3Ku/2 - 
3Vg U /2), and Ae gu = (Ae gu + 3X g /2 + 3X u /2). The inter- 
molecular transfer integrals t gg , t uu , and t gu are shown in Fig. 
[3] In terms of the basis given in Eq. (4.18l, this CI matrix can 
be block diagonalized as 



where £7 g ^ = 
3Vg U /2), £7< i 1) - 



(£7(°) + 

(£7(0) 
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H 2 . 



2-mol 



-Ae gu - 1 




gu 



-A5g U - 1 






-As 

M 1} 






gu + ^gu 
^1111 



(4.20) 



Since the MOs are frozen to the ones obtained in the 
isolated molecule calculation, the off-diagonal elements 

(V^i I #2-moi | ^ g ]l ) and (^l I ^2-moi | ^2 ) become non zero. 

Similar implementations can be performed for the 2- 
electron case. The expressions of each configuration on the 
molecular-based fragment MO basis and the corresponding 
energies in terms of the extended Hubbard Hamiltonian are 
shown in Table [TTJ We also analyze the off-diagonal compo- 
nents, and derive the CI matrix explicitly in terms of the pa- 
rameters of the extended Hubbard Hamiltonian. 



E. Results 

Now we can evaluate the model parameters by relating the 
CI matrix obtained by the MR-CI calculations given in Eq. 



(4.16) and the CI matrix expressed in terms of the model 
Hamiltonian. In the extended Hubbard Hamiltonian [Eq. 



(4.8)], there are 17 parameters, including the intra-molecular 
interactions. In the two-electron case, 22 diagonal terms and 
30 independent off-diagonal terms can be determined in the ab 
initio Hamiltonian ( |4.16| ). Interestingly, the intra-molecular 
interactions, U g , U u , U' , and Jh, are also obtained from the 
calculations for the two-molecule system. These values are 
consistent with the results in Table [J 

For [Au(tmdt)2], the inter-molecular interactions are listed 



in Table III From the data of transfer integrals, we see that 



TABLE III: Estimated parameters for the neighboring [Au(tmdt)2] 
molecules. All energies are in eV. 



direction 


[100] 


[111] 


[101] 


[211] 


[001] 


[011] 


^gg 


0.10 


0.10 


0.02 


0.01 


0.07 


0.01 


^uu 


0.12 


-0.19 


-0.05 


-0.02 


-0.11 


-0.02 


tg u 


0.00 


-0.14 


-0.03 


-0.01 


-0.09 


-0.02 


Vgg 


1.75 


1.31 


1.38 


0.97 


1.37 


1.13 




1.67 


1.39 


1.43 


1.00 


1.37 


1.14 




1.70 


1.35 


1.41 


0.99 


1.38 


1.13 


I 


0.43 


-0.29 


-0.24 


-0.11 


-0.13 


-0.09 




0.09 


-0.56 


-0.52 


-0.28 


-0.45 


-0.31 


x u 


0.05 


-0.58 


-0.52 


-0.29 


-0.43 


-0.31 


Ae g 


0.26 


-0.11 


-0.06 


-0.07 


0.03 


-0.01 


Ae u 


0.20 


-0.06 


0.00 


-0.05 


0.07 


0.02 


A£ g u 


-0.10 


0.21 


0.19 


0.04 


0.31 


0.12 



the system exhibits three-dimensional character, in agreement 
with DFT-based calculations PS The estimated parameters for 
As g , Ae u , and Ae gu are also listed in Table III and are com- 
parable to the energy difference (s® — e®) « 0.26 eV. For 
(TTM-TTP)l3, the inter-molecular interactions exhibit strong 
anisotropy, i.e., the inter-molecular interactions for the [001] 
molecule pair become largest compared with those for other 
pairs. 41 These features can be explained by noting that [001] 
is the stacking direction of TTM-TTP molecules and the two- 
molecule distance becomes much shorter for this face-to-face 
molecule pair. The parameters As g and Ae u for the ionic 
TTM-TTP+ molecules are large in contrast to those for the 
neutral [Au(tmdt)2] molecules. Since these values are almost 
identical Ae g ~ Ae u , we find that the g and u MOs are still 
quasi-degenerate in the crystal. 

It is known that the orbital exchange interaction term / ac- 
counts for dispersion interactions (i.e., van der Waals interac- 
tions) between pairs of molecules r 3 " 45 ! These correspond to 
instantaneous dipole-dipole interactions resulting from the lo- 
cal charge excitations gi — » ui and gn — » un [see Eq. (4.6)], 
with the excitation energy AE. The dispersion interaction can 
be evaluated as ~ —I 2 /AE, and the typical amplitudes of the 
dispersion interaction are in the range —0.1 to —0.01 eV. In 
the present systems, the dispersion interactions for nearest- 
neighboring molecules are —0.24 eV for [Au(tmdt)2]2 and 
-0.08 eV for [TTM-TTP - ] 2 . For the typical one-band sys- 
tem of the TTF molecules, we also evaluated the parameter 
/ = 0.19 eV and found that the dispersion interaction is very 
weak (~ —0.006 eV). Such an observation can be ascribed to 
the fact that the present systems are very polarizable due to 
extended MOs. 

In the way shown above, we determine all the possible 
model parameters uniquely for multi-orbital systems by tak- 
ing the advantage of the wavefunction-based ab initio calcula- 
tions. However, the resulting parameters shown here are bare 
values and the screening effects are not taken into account. A 
framework to project the effective model from large CI expan- 
sions has been developed where the electronic excitations in- 
volving valence and virtual MOs , i.e., s o-called dynamic cor- 
relations, are taken into accountP- 1 ! 12 * 42 ! With this scheme, one 
can access to accurate magnetic interactions, and simultane- 
ously, transfer integral t, Coulomb repulsion U, and the di- 
rect exchange coupling, with the inclusion of these screening 
effects. For atomic orbitals in metal complexes, the screen- 
ing effects to Coulomb repulsions are pronounced due to the 
localized orbital nature and also due to the presence of lig- 
ands, typically by a factor of 4 to 5P^ It has been shown 
that this scheme is relevant to the analysis of pure organic 



11 



materials JBED a nd the screening effects would be less impor- 
tant for larger molecules. 5 From the recent approaches based 
on DFT, where the inter-molecular screening effects in addi- 
tion to the above intra-molecular screening effects are taken 
into account, the bare magnitudes for the Coulomb repulsions 
are reduced to at most a quarter.33 



V. MODEL FOR CRYSTAL SYSTEM AND BAND 
STRUCTURE 

In this section, we consider the crystal system and evalu- 
ate the band structure. The symmetry-breaking term Ae gu 
[Eq. (|4.4|)] disappears in the periodic crystal system. The full 



Hamiltonian is given by H ( 



cryst 



rrkin 

cryst 



rrint 

cryst 



with 



kin 

cryst 



rrint 

cryst 



£(e g n g ,, + e u n uJ ) - £ £^'^(4,^ W + h - c " £ '""'"M.i^j,* + ^0 

3 (iJ) a a 

- ^^K,n b ,n c ]/ t . c . _ c t . i he) 



(5.1a) 



+ £ V^ n ^n sA n sJ + £ Vt- nb '^ ] n u ,in ud + £ T/K>™^] (n g>i n UJ + n uA n sJ ) 

(id) (hj) 



+ X) ^4 na,nb ' nc] [>* ( C L> C ^> + h.c.) - (cl ha c u ^ a + h.c.) n gJ 



(5.1b) 



where (i, j) denotes the combination of neighboring molecule 
pair and is assumed that i represents the molecule in a refer- 
ence position while j the translated molecule by the vector 
[n a , n b , n c ] shown in Fig. [4] 

Next we examine the band structure for crystals by ne- 
glecting the correlation effects. In order to determine the en- 
ergy levels, the crystal field effect arising from the Ae g and 
Ae u terms must be taken into account. We note that, for 
(TTM-TTP)l3, the potential due to the counterion I3 should 
also be included in order to examine this crystal field effect 
quantitatively. In addition, the Hartree corrections from Eq. 
( 5.1b| ) also contribute to the energy difference of MO level 
energies; however, these contributions would be small since 
the density operators in Hamiltonian are represented in the 
normal-ordered form [see Eq. ( |3.7| )]. In the present anal- 
ysis for [Au(tmdt)2], we assign the MO level energies as 



2A4 100] and e u * s" u + 2A^ 



100] 



where the fac- 



tor 2 reflects the coordination number, and we simply focus 



on the kinetic term [Eq. ( |5.1a| )]. The resulting band structure 
for the [Au(tmdt)2] crystal is shown in Fig. [5] The bandwidth 
obtained from the present analysis is overestimated in com- 
parison with the result of DFT-based calculation for the peri- 
odic system, 22 however, qualitative behavior of the band struc- 
ture is well reproduced. By taking advantage of the present 
scheme, we can elucidate the nature of the band structure, by 
setting the mixing term to zero, i.e., t gu = 0. Such a fictitious 



band structure is also shown in Fig. [5] We observe that the two 
bands overlap and mix in together by large mixing amplitude, 
supporting a multi-band system. 

The band structure for (TTM-TTP)l3 is also analyzed.^ 
From an extended Hiickel approach, 4 the overlap integrals be- 
tween the neighboring molecules along the stacking direction 
were given by S gg = -0.19 x 1(T 3 , S uu = -26.22 x 1(T 3 , 
and S gu = — 12.92 xlO -3 . The small transfer integral t gg and 
also the small overlap integral S gg result from the tilted align- 
ment of TTM-TTP molecules along the stacking direction in 
the crystal. We observe that the band built on the u MO is 
much wider than the one built on the g MO. The bandwidth 
for the u MO is - 1.2 eV while that for the g MO is ~ 0.3 eV. 
Finally, since the inter-orbital transfer integral tg° u 01 ^ (= —0.13 
eV) is relatively large compared with t gg , the g MO band is 
strongly modified, while the u band is not much affected near 
the Fermi energy. 



VI. FULL FRAGMENT DECOMPOSITION 

In this section, we derive the effective model based on full 
fragment MOs, which is the most fundamental model to an- 
alyze the "intra-molecular" degree of freedom. The relevant 
fragment MOs are simply the left and right part of MOs, 
and (/pr, where the center fragment is omitted due to its small 
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x r y r z s 



FIG. 5: Band structure of [Au(tmdt)2] obtained from a parametriza- 
tion based on our ab initio calculations, where F = (0, 0, 0), X = 

(tt, 0, 0), Y = (0, tt, 0), Z = (0, 0, tt), S = (tt, tt, 0), U = (tt, 0, tt), 
T = (0, 7r, 7r), and R = (-zr, tt, it). The dotted curves represent the 
energy dispersion setting £ gu = 0. 



weight. 29 We denote the corresponding annihilation operators 
by CL,j,cr and cr^. The operator correspondence between 
the original MOs and the fragment MOs basis is given by 



J2 J2 t[ 2 a,nb,UC \4l,i 7 * C RJ,<7 + C L,i,a C L,j,a + h.C.) 

EE^'" J " ,cI (4, i /R,i.+h.c) 

(ij) G 

E T r[n a ,nb,n c l 



[n a ,n b ,n c ] 



J2vt a ^ nc] n L ,in^ 

(iJ) 



J ' 



(6.3) 



where a =t (|) for a =| (|). The parameters Eq and to 
represent the energy level of the fragment MO and the inter- 
fragment transfer integral within the molecule. The parame- 
ter U represents the magnitude of the Coulomb repulsion be- 
tween electrons within the fragment MOs, and Vb, J, and X 
denote the inter-fragment Coulomb repulsion, exchange in- 
teraction, and bond-density interactions within the molecule, 
respectively. The inter-molecular transfer integrals U (i = 
1, 2, 3) and the inter-molecular Coulomb repulsions Vi (i = 
1, 2, 3) are depicted in Fig. [6] The density operators are given 
in the normal-ordered form n L/Rj> = (c[ /Rj> c L/R j (T - §) 
and the basis- set change leads to 



( — C L,J> + C R,j,cr) j 



(6.1a) 
(6.1b) 



Let us mention that (fL and <£r are not orthogonal, but the 
overlap integral is very smallP^ 

On the fragment MO basis, the model Hamiltonian for 
the crystal [Eq. \5A\ \ can be re-expressed as H cryst = 

rrintramol i rrintermol „ 7 ;+u 



H^™™ l + H^™* with 

rrintramol 

crystal 



3 

3 ° 
3 

VoY n ^3 n ^j 

3 

J Y Sh >j ' Sr >3 ^LvT C L\I C ^3,i C ^t + h - C 
3 

X ^(^LJ.o- + ™Rj>)( c Lj\* C Rj\* + h - C 0> ( 6 - 2 ) 



intermol 



rr 

crystal 

= -EE t i w,l ( c U'+ k ) 



(6.4a) 
(6.4b) 



= ^ g + E/u) + ^' + | Jh, (6.4c) 
Vo = ^(£/ g + CAi) + ^' - ^ Jk, (6.4d) 



^ = ^(t/g + ^-^-JjH, 



1 



X = -(-/7 g + /7 u ), 



(6.4e) 



(6.4f) 



for the intra-molecular parameters. This kind of transforma- 
tion was applied for simpler two-orbital systems. 47 Similarly, 
the inter-molecular interactions can be expressed as 



ti — 2 ( tgg H~ ^im H~ 2tg U ), 



t2 — -(^gg+^uu), 



^3 — 2^ ^gg + ^uu 2t gu ), 



(6.5a) 
(6.5b) 
(6.5c) 



Vi = \(V gg + V uu ) + ^V gu - I - X g - X u , (6.5d) 



^2 = J(^gg + Ku) + ^ gu +/, 



(6.5e) 
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FIG. 6: The full-fragment-MO picture and the definition of parameters for (TTM-TTP)I 3 (left) and [Au(tmdt) 2 ] (right). 



V 3 = -(^ gg + Ku) + -V gu -I + X g + X u , (6.5f) 

where the index [n a , n^, n c ] is suppressed. The evaluated pa- 
rameters are summarized in Table II VI It is worthwhile to note 
that the presence of the nontrivial orbital exchange interaction 
I and the bond-charge interaction X g and X u plays a crucial 
role to differentiate the Coulomb interactions, Vi, V2, and V3 
in the fragment-MO picture. 

The magnitudes of the Coulomb repulsion within the frag- 
ment, U ~ 5.30 eV for TTM-TTP+ and U ~ 5.26 eV for 
[Au(tmdt)2], are comparable with that for the TTF molecule 
Z7ttf — 6.2 eV. For the inter- fragment transfer integrals, not 
only the nearest-neighbor Coulomb interactions but the long- 
range interactions have large amplitudes. In order to verify 
the fragment decomposition, we examine the distance depen- 
dence of the inter- fragment interactions. We define the inter- 
fragment distance r by the averaged inverse distance between 
the TVs = 6 sulfur atoms in each fragment, by 

\ = w^h (6 - 6) 

where is the distance between the sulfur atoms. The 
r-dependences of the inter-fragment interaction V for 
[Au(tmdt) 2 ] and (TTM-TTP)I 3 are shown in Fig. [7] and in the 



supporting information. 41 The dotted curve denotes the bare 
Coulomb interactions V(r) = l/(47T£oO> wnere £ o is the per- 
mittivity of vacuum. Thus, we conclude that V follows well 
the Coulomb law, and our fragment decomposition is verified 
from this evaluation of inter- fragment Coulomb repulsion. In- 
cidentally, the fact that the intra-molecular Coulomb repulsion 
Vo also follows well the Coulomb repulsion supports that the 
center fragment can be neglected. 

From the data of inter- fragment transfer integrals, we find 
that the stacking TTM-TTP molecules ([001] direction) can be 
described as a two-leg ladder system, where the transfer inte- 
gral along leg direction is £2, while those along the rung direc- 
tion are to and t\. It is worth noting that the inter-molecular 
transfer integral t\ exceeds the intra-molecular transfer inte- 
gral to. This implies that, as a possible origin to the insulat- 
ing behavior at high temperature, electrons are localized on 
each t\ bond, a reminiscence of the "dimer-Mott" state in the 
quarter-filled ID systems ! 48 " 5Q I 

Here we discuss the low-temperature symmetry-broken 
states in (TTM-TTP)I 3 on the basis of the present model 
parameters. In (TTM-TTP)l3, the intra-molecular charge 
ordering has been proposed from the Raman and x-ray 
experiments E3EU However, the charge pattern and the origin 
of the spin- singlet behavior have not been clarified yet. From a 
simple strong-coupling analysis, we can examine the energies 
of possible ordered states as shown in Fig. [8] We observe that 
the lowest-energy state is the ICO state (£?ico = —1.47 eV 



TABLE IV: Estimated parameters for the inter-molecular interac- 
tions on the full-fragment MO basis for (TTM-TTP)l3 . The intra- 
molecular interactions are given by so = —12.07 eV, to = —0.17 
eV, U = 5.30 eV, V = 2.07 eV, J = 0.18 eV, and X = 0.05 eV. 
All energies are in eV. 



TABLE V: Estimated parameters for the inter-molecular interactions 
on the full-fragment MO basis for [Au(tmdt)2]. The intra-molecular 
interactions are given by £0 = — 5.70 eV, to = 0.19 eV, U = 5.26 
eV, Vo = 1.58 eV, J = 0.08 eV, and X = -0.01 eV. All energies 
are in eV. 
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FIG. 7: The Coulomb repulsions as a function of inverse inter- 
fragment distance r _1 for [Au(tmdt)2] . The results for (TTM- 
TTP)l3 are shown in the supproting information. 41 The filled, open 
and double symbols represent Vi, V2, and V3, respectively. The 
repulsion between the L and R fragments within the molecule is 
characterized by Vo . The Coulomb repulsion within the fragment is 
U ~ 5.26 eV. The dotted line represents the bare Coulomb repulsion 
a/r with a « 14.4 eV A. 



per molecule) in which the charge is disproportionated within 
each constitutive molecule. This ICO pattern is compatible 
with the qi = (0 . 0, 1/2 ) superstructure observed in the x- 
ray measurements .klE^ES From the evaluated model parame- 
ters, we expect that the super-exchange interaction along the 
t\ bond becomes largest and it would play a crucial role to 
induce the spin- singlet state. On the basis of this finding, we 
infer that the non-magnetic insulating behavior observed at 



low temperatures in (TTM-TTP)I^^can be attributed to the 
spin-singlet formation on the 4°°^ bond (shown by the dotted 
circle in Fig. [8]) with twofold periodicity along the stacking 
direction. Incidentally, we observe that the IAF state is almost 
degenerate with the ICO state, with the energy differences per 
molecule E\ KY - ^ico ~ 0.15 eV, E S bw ~ ^ico ~ 0.30 eV, 
and ^cdw — ^ico ~ 1.19 eV. Detailed analysis of possible 
symmetry-broken states described by the present Hamiltonian 
has been reported by using the mean-field approximation.^ 

From the optical measurement analysis, a nontrivial op- 
tical absorption band has been observed at « 5000 cm -1 , 
corresponding to a charge transfer band!^ It has been sug- 
gested that this band can only be observed under some elec- 
tric fields polarized perpendicular to the stacking direction 
E JL c. Furthermore its intensity is strongly enhanced at 
low temperatures. A possible scenario to explain this be- 
havior is the following. In the ICO state, a characteristic 
charge excitation perpendicular to the stacking direction can 
be described with two holes localized on two adjacent frag- 
ments (see dotted circle in Fig. [5]). If one considers a two- 
site two-electron system as a simplest model for this unit, the 
ground state is singlet with energy Eo = (U + v/ 001 ^)/2 — 
[(U - 1 V 1 [001] ) 2 /4 + 4(4 001] ) 2 ] 1 /2. while the first-excited state 
with energy E\ = V\ represents an optically-forbidden spin 
transition, the unique optically-allowed transition involves the 
second-excited state with energy E2 = U, representing a 
charge excitation. The corresponding excitation energy can be 
roughly estimated as 6800 cm -1 by using reduced Coulomb 
values (by a factor of 5). 58 Despite the localized character 
of this excitation description, our result is in relatively good 
agreement with experiments. 



ICO 




IAF 



CDW 



c 

o 

e 



SDW d>>*>>* 

1 8 8 / 8 8 T 




FIG. 8: Strong-coupling picture for several possible ordered states 
in (TTM-TTP)l3 shown on the fragment-MO basis. Each circle (ar- 
row) in the ICO and charge-density-wave (CDW) states represents 
one hole (spin), while that in the IAF and spin-density- wave (SDW) 
states represents half of hole (spin). In the ICO state, two holes sur- 
rounded by the dotted circle represent a spin- singlet pair. In the CDW 
state, two holes form a spin-singlet pair within a molecule. 



Similar analysis is performed for [Au(tmdt)2]. For 
[Au(tmdt)2], we also find that the inter-molecular transfer in- 



tegrals, t^ °\ ^ ±±±J , and £i UU±J exceeds the intra-molecular 
transfer integral to. These features are qualitatively consistent 
with those obtained by fitting the DFT-based calculation.^ 
In comparison to (TTM-TTP)l3, we found that the two 
fragments connect by the tj 111 ' bond form a strong dimer 



[nil 



,[001] 



(4 111] /4 m] 



7.0) in [Au(tmdt)2], while four fragments in- 



teract simultaneously in (TTM-TTP)l3 where t^ 01 '/^ 001 ' ~ 
1.5. Due to this feature, the ICO state becomes unfavorable 
in [Au(tmdt)2]. In addition, it has been pointed out that this 
system has a good nesting vector q = ( 1/2, , 0) and the 
IAF state with this wave vector is stabilized!^!] If we restrict 
ourselves to the ordering with wave vector q = (1/2,0,0), 
we observe from the strong-coupling analysis, that the ICO 
and CDW states are unstable with respect to the charge- 
uniform state in which each hole is localized on the t^ 11 ^ 
bond. We note that this charge-uniform state is compatible 
to the IAF state if we take into account the antiferromag- 
netic interactions. The energy differences per molecule are 
Eico-Ei AF « 0.07 eV and E C bw -#iaf ~ 0.14 eV. More 
elaborate calculations based on our effective model should be 
carried out to clarify the origin of a huge magnetic moment 
suggested by the nuclear magnetic resonance measurement.^ 
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VII. SUMMARY 

In the present paper, we have proposed a scheme to de- 
termine the parameters of multi-orbital extended Hubbard 
model from the ab initio MR-CI calculations. To the best 
of our knowledge, this is the first theoretical work which 
aims at evaluating model parameters for a multi-orbital sys- 
tem. We have applied this method explicitly to the charge- 
transfer molecular conductor (TTM-TTP)i3 and the single- 
component molecular conductor [Au(tmdt)2]. By taking ad- 
vantage of wavef unction-based calculations, the CI Hamilto- 
nian matrix for the target model space was constructed, and 
all the model parameters were uniquely determined so as to 
reproduce the different matrix elements. By examining the 
band structure, we have verified the multi-band nature of these 
systems, since the SOMO- and HOMO— 1 -based bands over- 
lap and these bands mix in together by the relatively large 
mixing amplitude. Furthermore, a full fragment decomposi- 
tion picture leading to a parameter hierarchization has been 
justified by the observation that the inter-molecular Coulomb 
repulsions as well as the intra-molecular interaction Vq follow 



well the Coulomb 1/r law. Our results strongly support that 
the ICO state experimentally-observed in (TTM-TTP)l3 must 
be described by a multi-orbital picture. 
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